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Abstract 

In this paper, we study the bifurcation of limit cycles in Lienard systems of the form ^ = 
y — F[x) , -j| = —x , where F(x) is an odd polynomial that contains, in general, several free 
parameters. By using a method introduced in a previous paper, we obtain a sequence of algebraic 
approximations to the bifurcation sets, in the parameter space. Each algebraic approximation 
represents an exact lower bound to the bifurcation set. This sequence seems to converge to the 
exact bifurcation set of the system. The method is non perturbative. It is not necessary to have a 
small or a large parameter in order to obtain these results. 
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The Lienard equation : 



d 2 x 



+ f(x)— + x = 



(1) 



appears very often within several branches of science, such as physics, chemistry, electronics, 
biology, etc (see [§, [| and references therein). 

This equation can be written as a two-dimensional dynamical system which reads as 
follows : 



number and location of limit cycles. 

In order to make progress with this problem, it is of fundamental importance to control 
the bifurcations of limit cycles that can take place when one or several parameters of the 
system are varied. The word bifurcation is used to describe any sudden change that occurs 
while parameters are being smoothly varied in any dynamical system. Connections with the 
theory of bifurcations penetrate all natural phenomena. The differential equations describing 
real physical systems always contain parameters whose exact values £1X6 ; cLS £L rule, unknown. 
If an equation modeling a physical system is structurally instable, that is if the behavior 
of its solutions may change qualitatively through arbitrary small changes in its right-hand 
side, then it is necessary to understand which bifurcations of its phase portrait may occur 
through changes of the parameters. 

In this respect, the most difficult bifurcation is the so-called saddle- node bifurcation of 
limit cycles : let us suppose that system (Q) depends on a parameter A : F = F(x,X). Let 
T be a non-hyperbolic limit cycle of (0) (see |J for a definition) , corresponding to the 
value A of the parameter A. System (||) undergoes a saddle-node bifurcation at A = A if 
for A = Ao + e and e positive and sufficiently small, the limit cycle T bifurcates into two 
hyperbolic limit cycles, one stable and the other instable. Moreover, for A = Ao — e, the 
limit cycle T disappears and there is no limit cycle in a small neighborhood of r . This 
bifurcation is particularly difficult to detect because for A = Ao — e there is no trace of 
it. Moreover, the value of Ao is not known in principle and it is not possible to employ a 
perturbative method with respect to e to study this type of bifurcation. 

In a previous paper 0, we have introduced a method for studying the number and 
location of limit cycles of (fj), for the case where F(x) is an odd polynomial of arbitrary 
degree. The method is as follows : we consider a function h n (x,y) given by : 




where F(x) = J* f(r)dr. 

The most difficult problem connected with the study of equation (Q) is the question of the 



h n (x, y) = y n + g n -x,n{x)y n + g n -2 



n 



(x)y n 2 + ... + gi >n (x)y + g ,n( x ) 



(3) 



where gj n (x), with j = 0, 1, n — 1, are functions of only x and n is an even integer. Then 

it is always possible to choose the functions gj t n{%) m such a way that : 

dh Oh 
h n {x, y) = {y - — - x- 



dx dy 

is a function only of the variable x (see also M). Then we have : 

h n (x,y) = R n (x) (4) 

The functions gj, n {x) and R n {x) determined in this way are polynomials. As explained in 
||, if for a given value of n, the polynomial R n {x) has no real roots of odd multiplicity, then 
the system has no limit cycle. 

We want to show in this paper that the method presented in || enables us to determine 
algebraic approximations to the bifurcation sets of limit cycles for the Lienard equation. 
These bifurcation sets can be determined analytically only when the system has a small 
parameter or a large one (perturbative regime). In the intermediate case (non-perturbative 
regime), no method is known for determining, in an analytic way, the bifurcation set. We 
shall show here that our method gives a sequence of algebraic lower bounds to the bifurcation 
sets. Moreover, this sequence seems to converge to the exact bifurcation set. The method 
can be applied to any system (Q) where F(x) is an odd polynomial. 

As an example, we will consider a Rychkov system : 

F(x) = ci2X 5 + aix 3 + aox (5) 

with a 2 7^ 0. 

We can take one of the parameters equal to one without loss of generality. Several authors 
have studied this system with F(x) written as F(x) = e(x 5 — /ix 3 + x). Rychkov has shown 
in H that this system can have at most two limit cycles and actually has exactly two limit 
cycles when e > and fi > 2.5. Rychkov's results have been improved by Alsholm 0, who 
lowered the bound of \i to 2.3178 and by Odani [0, who obtained an even smaller value \/5. 
By a scaling of the variables x and y, system (0), with F(x) given by (|5D, can be written in 
a more simple form, as follows : 

F(x) = x 5 - fix 3 + 5x (6) 

Since there are two parameters, the bifurcation set is given by a curve in the parameter 
plane (/i,5). O ur aim, here, is to obtain information about the bifurcation diagram of the 
system in this plane. 



• For 5 < 0, thanks to Lienard theorem (see 0), we know that the system has exactly 
one limit cycle for arbitrary values of fi. 

• For 5 > and \i < 0, the Bendixon criterium (see ||) enables us to conclude that the 
system has no limit cycle (the divergence of the vector field, given by —F'(x), has a 
constant sign for all x). 

• For 5 > and fi > 0, the system can have two or zero limit cycles , according to 
Rychkov's results. In this region of the parameter space there exists a bifurcation 
curve B(fi, 5) = 0. In the region B(fi, 5) > 0, the system has exactly two limit cycles 
and for B(fi, 5) < 0, the system has no limit cycle. On the curve B(fj,, 5) = the system 
undergoes a saddle-node bifurcation : there is a unique non-hyperbolic (double) limit 
cycle. 

Obviously, the function B(fi, 5) is not known and no analytical method for obtaining this 
function for arbitrary fx and S exists. We shall obtain a sequence of algebraic approximations 
to the function B(/i, 5). 

For a given even value of n, let us consider the corresponding polynomial R n (x). The 
polynomials R n (x) described above, can have, for system @, one, two or zero positive simple 
roots, depending on the values of fx and 5. At least it is the behavior observed for the values 
of n we considered. 

• For 5 < and V// , the polynomials R n (x) have one simple positive root. 

• For \i > and 6 > , the first quadrant is divided in two regions by a curve B n ((x, 5) = 
0. In the region B n (/j,, 5) > 0, the polynomial R n (x) has two positive simple roots while 
in the region B n (fi,5) < it has no positive root. On the curve B n (fj,,5) = 0, R n (x) 
has a double positive root. 

• For \i < and 6 > 0, the polynomials R n {x) have no real root other than the even- 
multiplicity root in x = 0. 

It is clear (see 0) that for fx > and 5 > lying in the region B n (fi, 5) < 0, the system (H) 
with F(x) given by (^ has no limit cycle. Hence, it is evident that the curve B n (fi,5) = 
represents an exact lower bound to the bifurcation curve B(fi, 6) = : the curves B n (fi, 5) = 
are contained in the region B(fi, 8) < for all even values of n. 

The functions B n (fi, 5) are algebraic and can be determined from the conditions : 

dR 

R n {x) = and -^{x) = (7) 



These two algebraic equations determine the double root of the polynomial R n (x) and give 
a relation between \i and 5 which we write B n (fi,5) = 0. For n = 2, we find B 2 (fi,5) = 
5(fi 2 —4(5). For n = 4, B^/i, 5) is a 12 th degree polynomial. The degree of -£>„(/i, 5) increases 
rapidly with n. We have calculated the functions _B n (/i,<5) for even values of n between 2 
and 14. The behavior of the curves B n (fi, 5) = 0, as well as the numerical bifurcation curve 
(calculated from a numerical integration of the system), are shown in fig. ([!]). We see that 
each curve B n (/i, 5) = is contained in the region B n ^ 2 (n, 5) > 0. The complete bifurcation 
diagram is given in fig.(^). There are three regions : 

Region I : There is no limit cycle. All the curves B n (n, 5) = lie in the part /i > of 
this region. 

Region II : There are two limit cycles. All the polynomials R n {x, /i, S) have two positive 
simple roots. 

Region III : Lienard theorem shows that there is one limit cycle. All the polynomials 
R n (x,fi,5) have one positive simple root. 

We would like to emphasize that the shape of the bifurcation curve B(fi, S) = is already 
given by the curve B2(fi,5) = 5(/i 2 — 45) = 0, which is constructed only with the function 
F{x) ! The Hopf bifurcation happens when 5 = and the saddle-node bifurcation occurs 
near fi 2 = 45. 

We shall now make use of this bifurcation curve for the following system : 

x = y-(x 5 — (1 + A)x 3 + Xx) 

y = -x (8) 

This example has been studied by Lloyd || and Perko |J. System ([8]) is a particular case 
of system with : 

fi = -^-(1 + A) and 5 = X 
o 

In order to know what are the bifurcations of this system when the A parameter is varied 
from — oo to +oo, we must plot the line : 

A* =^(1 + 5) (9) 

in the bifurcation diagram of system (|j), with the exact (but unknown) bifurcation curve 
B(n,5) = replaced by one of the algebraic approximations B n (fi,5) = 0. As A is varied 
from — oo to +oo, the system moves along the line (|§D from left to right in fig.(^). It is easy 



to see that when A is negative (the portion of the line is in region III), there is one limit cycle. 
Then, when the line crosses the fx axis (that is when A changes sign) the system undergoes a 
Hopf bifurcation : a small limit cycle is created around the origin of the phase plane. There 
are now two limit cycles, the system is in region II. But when A is further increased, the line 
(0) crosses the bifurcation curve B(fi, 6) = : the two limit cycles collapse in a saddle-node 
bifurcation and there is no limit cycle , the system is in region I. If we continue to increase 
A, we see the line @ crossing the curve B(fx, 5) = again : two limit cycles appear in a 
saddle-node bifurcation, the system enters region II again. We can see that this is the last 
bifurcation we can create because the line, when A is further increased, does not cross the 
bifurcation set any more and stays in region II. From the intersections between the line (^) 
and the curves B n (fx,5) = 0, we obtain algebraic approximations to the bifurcation values 
of the parameter A. 

There is another way to see the different bifurcations of (|j). We have claimed in || 
that if n is large enough, the number of positive roots of odd multiplicity of R n (x) gives the 
number of limit cycles of the system. In the case of system (H), R n = R n (x, A), so when 
A is varied, the number of roots of R n changes. Hence, for a given value of A, the number 
of limit cycles of (§) can be obtained by counting the number of intersections between the 
curve R n (x, A) = and the line A = cte in figure ([|). 

In ||, we claim that the value of the root of R n {x) gives an approximation to the 
maximum value of x on the limit cycle (which we call the amplitude). In fig.(^) we have 
plotted R n (x, A) = for n = 2 , n = 6 and n = 10. So we can see the amplitude of the limit 
cycles with respect to A. We see the Hopf bifurcation when A crosses the x-axis upward and 
we see the two saddle-node bifurcations when R n (x,X) loses its two positive roots. Once 
again, as R^fe, A) = — 2xF(x, A), an approximation to the bifurcation amplitude-diagram is 
given by the curve F(x, A) = 0, which can be written as : 



A = x 



2 3x 2 - 
V5x 2 - 3 



Let's consider another example : 



x 



y - (x{x 2 - a 2 ){x 2 - 2 2 ){x 2 - 5 2 



)) 



y = -x 



(10) 



We want to study the bifurcations of (|T0|), when a is varied from — > oo. As we already 
twice noticed, the qualitative bifurcation amplitude-diagram seems to be given by F(x) = 0. 



Here, the plot of F(x, a) = seems to announce the presence of a transcritical bifurcation 
(see p| for a definition) near the values a = 2 and a = 5 and a Hopf bifurcation for a = 
(see fig.@). If we plot R 4 (x, a) = 0, we still see the Hopf bifurcation, but the supposed- 
transcritical bifurcations are indeed saddle-node ones (see fig.@) : we see that the system 
can have one or three limit cycles. When a is far from the values a = 2 or a = 5, there 
are three limit cycles, but when a is near the values a = 2 or a = 5, there is only one limit 
cycle. So, in this example, the equation F(x, a) = does not give the right qualitative 
amplitude-bifurcation diagram. We must plot the curve i? 4 (x, a) = in order to obtain the 
good qualitative shape of it. 

In summary, we have introduced a method that gives a sequence of algebraic approxi- 
mations to the bifurcation sets of limit cycles for the Lienard equation (§). These algebraic 
approximations are exact lower bounds to the exact bifurcation sets of the system and seem 
to converge to it in a monotonous way. The fundamental aspect of this method is that it is 
not perturbative in nature. It is not necessary to have a small or a large parameter in order 
to apply it. 
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Figure 1: The algebraic curves B n (fi,5) with n = 2,6, 14 (continuous lines) and the bifur- 
cation curve B(/i, S) (dashed line) which is calculated from numerical integrations of system 




Figure 2: The complete bifurcation diagram of system (0) with F(x) given by @. On the 
line 5 = the system undergoes a Hopf bifurcation. On the curve B(/j,, 5) = we have a 
saddle-node bifurcation of limit cycles. The system has no limit cycle in region I, two limit 
cycles in region II and one limit cycle in region III. 
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Figure 3: The line fj, — ^(1 + 5) and the bifurcation curve -£>(/i, 5) — of system (^j. The 
intersections between the two curves give the bifurcation values of A for system (|]). 
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Figure 4: The algebraic approximations to the amplitude-bifurcation diagram for system 
(^|) for the cases n = 2 (bold), n = 6 (dash) and n = 10 (continuous). We see the number 
and the amplitudes of the limit cycles by considering the intersections between one of these 
curves and a line A = cte. The results are improved with increasing values of n. 
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Figure 5: The curve F(x,a) = for system (pT0|) . We see the number and the amplitude 
of the limit cycles by considering the intersections between this curve and a line a = cte. 
There seems to be transcritical bifurcations for a = 2 and a = 5 at this order. 




Figure 6: The curve R/±(x,a) = for system fllOl). This curve is an approximation to the 
amplitude-bifurcation diagram. We see the number of limit cycles by counting the number 
of intersections between this curve and a line a = cte. The system only presents saddle-node 
bifurcations. The qualitative behaviour of the curves R n (x, a) = 0, with n > 4, is the same. 
No further qualitative changes occur for greater values of n. 



